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A path integral formulation is used to study the fluctuation-induced interactions between mani- 
folds of arbitrary shape at large separations. It is shown that the form of the interactions crucially 
depends on the choice of the boundary condition. In particular, whether or not the Casimir in- 
teraction is pairwise additive is shown to depend on whether the "metallic" boundary condition 
corresponds to a "grounded" or an "isolated" manifold. 



I. INTRODUCTION AND SUMMARY 



External objects that are immersed in a fluctuating 
medium, and modify the fluctuations in their vicinity, 
experience induced interactions with one another pj-p| . 
These interactions are most often independent of the 
structural details, and are in turn highly sensitive to the 
geometry of the objects and their mutual arrangements 
while immersed in the medium. 

The strong dependence of these interactions on the 
shape of the objects raise the issue of pairwise additivity: 
Is it possible to express the fluctuation-induced interac- 
tion between two extended bodies as the sum of a pair 
potential, or the interaction between several bodies as 
the sum of two-body interactions? 

It is well known that a pairwise summation of the van 
der Waals interaction gives the correct power law for the 
Casimir energy j^]. Let us take a pair potential of the 
form — A/r n , with n = 6 for the thermal case and n = 7 
for the quantum case and A being a constant to be 
determined. If one tries to fix the coefficient by summing 
the pair potential over two bodies and equating the result 
to the expression for the Casimir interaction between the 
bodies, one finds out that a different coefficient is needed 
for every geometry. 

To understand this, one should note that the van der 
Waals interaction is due to dipolar fluctuations. When 
two extended bodies are at a close separation, one can 
show that the fluctuations of all the multipoles in fact 
contribute comparably to the Casimir energy, and thus 
summation of the contribution due to the dipolar fluc- 
tuations cannot by itself account for the interaction || . 
When the bodies are at large separations (larger than 
their typical sizes), the contribution due to higher multi- 
poles are in fact systematically weaker. However, there is 
still a discrepancy between the sum of a (van der Waals) 
pair potential, and the contribution of the dipolar fluc- 
tuations to the Casimir energy. In the spirit of a (second 
order) perturbation theory, the correct way of calculat- 
ing the dipolar Casimir energy is to consider the pairwise 
sum of the dipole-dipole interactions over the two bod- 



ies, and then square it and take the average. This is 
clearly in contrast with the pairwise summation of the 
van der Waals interaction, which corresponds to taking 
the square of the local dipolar fluctuations and averaging, 
and then summing over the two bodies. 

The same picture can help us answer the second ques- 
tion. Many-body interactions can be expressed as the 
sum of many-body interactions of the multipoles of dif- 
ferent bodies in the medium. When extended bodies are 
at close separations, and all the multipoles have compa- 
rable contributions, many-body interactions of nontrivial 
forms result (7|-[|. On the other hand, for bodies at large 
separations the leading order contribution comes from 
the sum of two-body interactions of the lowest nonvan- 
ishing multipole ]l0| . 

In this article, we study the issue of orientational 
pairwise additivity pCfl , which is to determine whether 
the orientational dependence of the interactions could 
be obtained from the summation of a pair poten- 
tial. A path integral formulation is used to study the 
fluctuation-induced interactions between manifolds of ar- 
bitrary shape at large separations, in the context of a 
multipole expansion. It is shown that the form of the in- 
teraction crucially depends on whether the manifolds are 
grounded or isolated in an electrostatic analogy. In the 
grounded case, the manifolds are connected to a charge 
reservoir to maintain a constant potential, and thus the 
leading fluctuations are monopolar. Isolated manifolds, 
however, are constrained to have fixed overall charges, 
and can only undergo dipolar fluctuations. The leading 
interaction between grounded manifolds is found to be of 
the form (monopole — monopole) 2 , and is independent of 
their shapes and orientations. The leading shape depen- 
dent term comes from the (monopole — dipole) 2 term, 
which gives rise to orientational dependencies that are 
pairwise additive. The interaction between isolated man- 
ifolds, however, is dominated by the (dipole — dipole) 2 
term to the leading order, which is not pairwise additive. 

The rest of the paper is organized as follows. In Sec. 
IJ, the path integral formulation is developed and gen- 
eral expressions are derived for the fluctuation-induced 
interactions for different types of boundary conditions. 
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In Sec. [II, the interactions are examined for the specific 



examples of symmetric objects such as spheres, and also 
highly asymmetric objects such as rods and disks, where 
the above features can be manifestly understood. Criti- 
cal fluids are examined in Sec. fw[ as a special case, and 
a conclusion follows in Sec. [v|. 



II. PATH INTEGRAL FORMULATION 

Consider a d-dimensional medium, in which a field <p 
is undergoing scale-free (massless) thermal fluctuations, 
described by the Hamiltonian 



H[<f>] = 



K 



(1) 



The field could represent a component of the electromag- 
netic field in a dielectric medium || , the electrostatic po- 
tential in charged fluids at very low salt concentrations 
te| , an order parameter field for a critical binary mixture 



or a magnetic system |TI]] , or a massless Goldstone mode 
arising from a continuous symmetry breaking |Q. For 
simplicity, in what follows we are going to think in the 
context of charged fluids, and later on comment on the 
specific effects in the case of critical fluids. 

Let us assume that there are n manifolds immersed in 
the medium, denoted by M a (a = 1, • • • , n), which mod- 
ify the fluctuations. In the electrostatic context, one can 
view each manifold as a conductor that requires a con- 
stant value for the potential field in the whole volume 
that it encloses. A restricted partition function, which 
requires a value of (j) a for the potential field on the crth 
manifold, can then be written as 



~ n 

Z[<f> a ] = J 2ty(x) l[S{<P\ 



j\/„ 



-hi 



(2) 



Following Ref. J/J , the functional delta functions can next 
be represented by introducing the Lagrange multiplier 
fields p«(x), as 



Z[<f> a \ = J 2ty(x) JJJ^ V Pa (x) cxpj-y J d d x(V0) 2 +i^ J d d xp Q (x)[0(x)-0 Q ]| 

= Z x f[ f Pp Q (x) exp J I d d xd rf x'p a (x)G(x - x')^(x') -iV> a / d d xp a (x) \ 



(3) 



in which 

G(x - x') = 

with Sri = 



(-v 2 ) 



2^-1 
x.x' 



1 



5 d (d-2)|x-x' 



/Id— 2 : 



(4) 



T(d/2) 



(the surface area of the d- 
dimensional sphere), Zq is the free partition function, 
and J M T>p a {yi) implies a functional integration only in 
the region enclosed by M a . (In other words, the Lagrange 
multiplier field /o a (x) is nonzero only within the volume 
of M a .) Note that one should view the p a (x) fields as 
fluctuating charge density fields, and Eq.(||) as the par- 
tition function of a set of interacting Coulomb plasmas, 
in the electrostatic context |l2| . 

The fluctuation-induced interactions between the con- 
ductors can now be inferred from the above partition 
function. However, it is important to specify the bound- 
ary conditions for the conductors. One possibility is that 
the conductors are grounded, that is to say they are main- 
tained at a constant fixed potential by being in contact 
with a large reservoir of charges; a so-called "ground." 



d d xd d x'- 



In this case, the free energy of the system is obtained as 



F r 



-k B T\nZ[(j) a = 0]. 



(5) 



The other possibility is that the conductors are made 
isolated, and maintain constant amounts of net charges, 
which we assume to be zero. In this case, the potential 
field at the conductors can take any value to help main- 
tain the neutrality, and thus the free energy is obtained 
as 



" + 00 



F is = -k B Tln[ / J\d<p a Z[<p o 



(6) 



To further proceed, we focus on the situation in which 
the manifolds are far from each other, namely, they are 
at separations much larger than their typical sizes. In 
this case, we can perform a multipole expansion for the 
charge density distribution. For example, the Coulomb 
interaction between the ath and the /3th conductors can 
be written as (a ^ 0) 



S d (d-2)|x-x' + R Q/3 | d - 2 ' 

QaQf} I QpPa ■ R-a/3 — QqP/3 ' R-a/3 



P Q • P/3 - dP n K^bPb R 



S d (d-2)R d - p 2 



Q T)d—1 



SdR d aP 



(7) 
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in which R a p is the distance between the two conductors, 
and the multipoles are defined as 



Q a = J rf d x/? Q (x) = p a (k) 
P a ,i = J d d xxip a (-x.) = ^--^-p a (k) 



k=0 



with 



p a (k) = J d d xp Q (x)e* kx . 



(8) 
(9) 

(10) 



We also need to make a similar multipole expansion 
for the self energy terms at each manifold 



d d xd d x' 

d d k 1 

(2ir) d ¥ 



p a (x.)p a (x.') 



5 d (d-2)|x-x 
p a (k)p a (-k). 



/Id— 2 ' 



(11) 



We can introduce the multipoles, using the Taylor ex- 
pansion of the charge density in Fourier space 



- /, ,s - / n s . 1 9p a (k) 
^( k )=A.(0)+ ? -g j5 - 



k=0 



(12) 



The above expansion can be formally viewed as an ex- 
pansion in powers of kL a , where L a is a typical size of 
the manifold. The expansion is thus convergent only for 
sufficiently small values of k, corresponding to length 
scales larger than the size of the manifolds. Since the 
self energy integral in Eq. (fTD) involves contributions from 



\{dQ c 



dV r 



higher wavevectors, a multipole expansion for the self en- 
ergy will be divergent. However, the expansion in Eq.jT^) 
indicates that all the information concerning the first few 
multipoles of the charge distribution is already contained 
in the low k behavior of the function p a (k) . Since we are 
only interested in the dependence of the partition func- 
tion Eq.(|3|) on the distances R a p> all we need to know 
about the self energy is its dependence on the first few 
multipoles, which is in fact well behaved. 

Let us denote the domain of convergence for the expan- 
sion m Eq.@ in k -space by T) a . This domain contains 
the origin, and its shape is determined by the geometry 
of the conductor. Loosely speaking, its size in each direc- 
tion is set by the inverse of the size of the conductor in 
that direction. Now we can restrict the fc-integral in the 
self energy only to this domain, and neglect the contribu- 
tion from the outside of T> a , because all the dependence 
on the first few multipoles is included in the domain T> a . 
Wc thus have 



d d k 1 



p a (k)p a (-k) 



) Va (2v) d k 2 

^)aO a ~t~ *Ya,ij P(x,iP(X,j 



in which 



7a 



7a 



d d k 1 

d k k j k j 
n a (2ir) d ~P~ 



(13) 

(14) 
(15) 



and so forth. 

Putting all the pieces together, the i? Q( g-dependent 
part of the partition function can be written as 



x exp ■ 



— E 

2K ^ 



(16) 



—y 

2K ^ 



QaQ/3 

S d (d-2)R d 



RaP,i{QpPa,i ~ Q a Pf3,i) ($ij ~ dR a p,iR a f3j)P a ,iP f 



a/3 



SdR 



d-1 

a/3 



S d Rt 



a/3 



Note that we have neglected a Jacobian in changing the 
measure of integration. However, since the transforma- 
tion from the charge density distribution to the multipole 
description is linear, one can show that the Jacobian is 



just an uninteresting constant. 

Finally, using Eqs.([|) and (|l6|), the interaction free en- 
ergy for grounded manifolds can be obtained as 



F„ = 



k B T 



E 

a#/3 L 



7a l 7/9 1 



(7a 



L 7/3 }j + 7 Q ,y7/3 l )R a p,iR a f} t . 



2 1? 



2(d-2) 



R 



2(d-l) 



0(1/ R 2d ) 



(17) 



The first term in Eq.(0) is a squared monopolc-monopolc interaction, and is independent of the relative orientations of 
the conductors in space. The second term, on the other hand, has the form of a squared monopole-dipole interaction, 
and does depend on the orientations through an effective dipole-dipole interaction, which is pairwise additive. 
Similarly, for isolated manifolds, Eqs.mfi and (Jlq) yield the interaction as 



k B T 



E 

a#/3 



7a,ifc7/3,j7 
K a0 



(Sij — dR a /3 } iR a p,j)(Ski — dR a f3^kRap,i) + 0(1/ R 



2d+2\ 



(18) 
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Note that the leading term in Eq.(18) is a squared dipole- 
dipole interaction, and thus it is not orientationally pair- 
wise additive. 



III. APPLICATION TO SPECIFIC GEOMETRIES 

The multipole expansion allowed us to calculate the 
general forms of the fluctuation-induced interactions be- 
tween manifolds of arbitrary shape and with arbitrary 
orientations with respect to one another, for the two cases 
of isolated and grounded boundary conditions. All the 
specific informations about the shapes and the orienta- 
tions of the manifolds are encoded in the 7-tensors de- 
fined above. These informations are in fact of three kinds: 
(i) the overall magnitude of the tensors which are set by 
the typical sizes of the manifolds, (ii) the orientational 
dependencies which make up the tensorial structure, and 
are dictated by the structure of the symmetry axes or 
"the principal axes" of the manifolds, and (iii) overall nu- 
merical prefactors of order unity. In this section, we try 
to use symmetry arguments to determine the 7-tensors 
for some simple geometries within the numerical prefac- 
tors, without actually specifying the exact shape of the 
integration domain T>. The final piece of information, 
which is the numerical prefactor, appears to be very sen- 
sitive to the exact geometry of the manifold (and thus 
to that of V), and can in general be calculated using the 
techniques developed in Ref. ECl). 



A. Two Spheres 



The 7-tensors for a sphere of radius L can be easily esti- 
mated using symmetry: 7 S ~ J^ L k d ~ x dkjk 2 ~ 1/L d ~ 2 



and jsAj ~ Sij J Q 



l/L 



,d-l 



dk 



5ij/L d . Using Eqs.rtl 



and another sphere of radius I/ 2 which is at a distance 
R, reads 



rd-2 t d-2 
KBJ X R 2(d-2) ' 



for the grounded case, and 



F 



,sph 



L d Li 



(19) 



(20) 



for the isolated case, with no orientational dependence 
due to symmetry. 



B. Two Rods 

The calculation of the 7-tensors for a rod of length L 
and thickness a is more tricky. Using the cylindrical sym- 
metry one obtains: "f r ~ f^ L dk z J^ a dk±k d l ^ 2 /(k 2 + 
k\) ~ 1/L d - 2 for d < 3, and - l/(La d - 3 ) for 
d > 3, where z-axis is parallel to the director of 
the cylinder, and _!_ denotes the remaining directions 
that are perpendicular to it. The second rank tensor 
"fr t ij is diagonal with only two independent components: 

~ ti /L dk z fc /a dk ± k<t~ 2 k 2 z /{k 2 z + k\) ~ 1/L d for 
3, and ~ l/(L 3 a d ~ 3 ) for d > 3, and 7 r .±j_ ~ 



7r,zz 

d 



< 
,i/z 



and rtlq), the interaction between a sphere of radius L±, 



So dkz So dk±ki/(k 2 + k 2 ) ~ l/(La d - 1 ). If the unit 
vector d denotes the director of the rod, the inverse sec- 
ond rank 7-tensor which appears in the expression for the 
interaction can be written as 7"^- ~ L d didj for d < 3, in 
the limit of small thickness. Note that in this limit, the 
inverse 7-tensors are vanishing for d > 3, and thus rods 
do not interact in these high dimensions. 

Using Eqs.(p7|) and (18), the orientation dependent 
part of the interaction between two rods of lengths L\ 
and L2, and directors d\ and d2, which are a distance R 
apart, reads (d < 3) 



prod 
gr 



T d-1 T d-1 

_fc B T x x 



R 2{d-l 



Li 
L 2 



(di • R12) 2 



(d 2 ■ R12) 2 



(21) 



for the grounded case, which is pairwise additive, and 



L d L d 

« x x 



d 1 -d 2 -d(d 1 -R 12 )(d 2 -R 12 ; 



(22) 



for the isolated case, which has a squared dipolar form 
and is not pairwise additive. 

C. Two Disks 

The 7-tensors for a disk of radius L and thickness a can 
be similarly calculated within numerical prefactors using 
symmetry. The zeroth rank tensor can be calculated as: 



Id - Jp 1/Q dk * Io 1/L dk ± k d ± - 2 /(k 2 z + k\) ~ l/L d ~ 2 , where 
z-axis is normal to the disk, and _L denotes the remaining 
directions in the subspace of the disk. The second rank 
tensor jd.ij is diagonal with only two independent com- 
ponents: 7 dj22 - J 1/Q dk z J 1/L dk ± k d ^ 2 k 2 /(k 2 + k\) ~ 

l/iaL*- 1 ), and ld>±± ~ J^ /a dk z f^ /L dk ± k d J \k 2 z + 
kj_) ~ 1/L d . If we denote the unit vector perpendicu- 
lar to the disk by n, the inverse second rank 7-tensor 
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which appears in the expression for the interaction can part of the interaction between a disk of radius L\ and 

be written as 7j",y ~ L d (Sij —hihj), in the limit of small normal vector ni, and another one with radius L 2 and 

thickness. normal vector n 2 that is a distance R apart, reads 
Using Eqs.(^7|) and (18), the orientation dependent 



jd-l rd-1 

£ isk ~ -k B T x L, 2 s x 



Li 

L-2 



l-(n!-Ri 2 ) 5 



L2 
Li 



1 - (n 2 • Ri 2 ) 2 



for the grounded case, which is pairwise additive, and 

d 2 - d - 2 + (ni ■ n 2 ) 2 + (2d - d 2 )(rn ■ R 12 ) 2 + (2d - d 2 )(h 2 • R 12 ) 2 
-2d(ni • R X2 )(n 2 • R u )(ni • n 2 ) + d 2 (ru • Ri 2 ) 2 (n 2 • R 12 ) 21 



F: 



disk 



k B T x 



T d t d 
R 2d 



(23) 



(24) 



for the isolated case, which has a squared dipolar form 
and is not pairwise additive. 



IV. CRITICAL FLUIDS 

As mentioned above, interactions could be induced be- 
tween objects that modify thermal fluctuations of an or- 
der parameter field for a critical binary mixture or a 
magnetic system ]Tlj . In this case, two kinds of bound- 
ary conditions are usually considered: (i) the ordinary 
boundary condition that suppresses the order parameter 
at the boundary, and thus does not break its symme- 
try, and (ii) the symmetry breaking boundary condition, 
which sets a nonvanishing value for the order parameter 
at the boundary. Note that the ordinary boundary con- 
dition is the same as the grounded boundary condition 
in the electrostatic terminology. 

The fluctuations in a critical fluid are characterized by 
the universality class of the system. For the case when 
the fluid can be described by a Gaussian Hamiltonian as 
given in Eq.(|l|), all of the above results for the grounded 
manifolds hold for the case of ordinary boundary condi- 
tion. The interaction between manifolds with symmetry 
breaking boundary conditions, where the value of the or- 
der parameter is set to $ Q on the ath manifold, is calcu- 
lated as 

F Bh = -k B T\nZ[(f> a =$J, (25) 
where Z[<f> a ] is given by Eq.(|i~6|). One obtains [O 



T^Gauss 



knT 



E 



S d {d-2)R d -* 



0(1/ R 



2d-4\ 



Eq. (|l7|) . This interaction is orientationally pairwise ad- 
ditive. 

For a nontrivial universality class, one should make use 
of more complicated Hamiltonians with nonlinear terms. 
It is then possible to calculate the Casimir energy ex- 
pressions using field theoretical techniques Q. The in- 
teraction between two spheres in an arbitrary critical sys- 
tem has in fact been calculated exactly in Ref. [[l4| us- 
ing conformal-invariance methods. The interaction for 
the case of symmetry breaking boundary conditions (on 
both spheres) is obtained as \/R d ~ 2+rt [ fTs] , while for 

the case of ordinary boundary conditions it is found as 
l/R2(d-i/v)^ where 7] and 

v are critical exponents of the 
system [0. One can easily check that for the case of 
Gaussian universality class, where 77 = and v = 1/2, 
they coincide with the results of Eq. ( |2^ ) and ( fl7| ) . 

It is interesting to note that the power law for the sym- 
metry breaking case is given by the two-point correlation 
function of the field (the spin-spin correlation in mag- 
netic terminology), while the one for the ordinary case is 
given by the four-point correlation function (the energy- 
energy correlation) [jl4|. Guided by this, one can think 
of an effective Gaussian Hamiltonian of the form 



K 
~2 



d d q 2 _ 
(2n) dq 



(27) 



(26) 



which yields a correct form for the two-point function, 
and calculate the fluctuation-induced interactions using 
the above methods 0. However, although it yields a 
correct result for the symmetry breaking case (almost by 
construction), it gives a corresponding form for the or- 
dinary case as l/_R 2 ( d-2 +'?) which is not correct. The 
reason is that the above effective Gaussian Hamiltonian 
does not give a correct four-point correlation function. 
However, it can be constructed to do so by using q 1 /" 
instead of q 2 ^ in Eq.d27 



It is important to note that this interaction is indepen- 
dent of the orientations of the manifolds, and that the 
leading order orientation dependent term for the sym- 
metry breaking boundary condition is the same as the 
case of ordinary boundary condition, and is given as in 



V. CONCLUSION 

The analysis that is presented here is aimed at em- 
phasizing the crucial role of the type of boundary con- 
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ditions on fluctuation-induced interactions. Unlike the 
case of extended objects at close separations, where dif- 
ferent types of boundary conditions all lead to the same 
form of interaction, we found that for objects at large 
separations the type of boundary conditions determine 
the very form of the interaction, and whether or not it is 
pairwise additive. 
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